Driven Dynamics of Periodic Elastic Media in Disorder 
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We analyze the large-scale dynamics of vortex lattices and charge density waves driven in a dis- 
ordered potential. Using a perturbative coarse-graining procedure we present an explicit derivation 
of non-equilibrium terms in the renormalized equation of motion, in particular Kardar-Parisi-Zhang 
non-linearities and dynamic strain terms. We demonstrate the absence of glassy features like di- 
verging linear friction coefficients and transverse critical currents in the drifting state. We discuss 
the structure of the dynamical phase diagram containing different elastic phases very small and very 
large drive and plastic phases at intermediate velocity. 
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I. INTRODUCTION 

Periodic structures driven through a random environ- 
ment have become a paradigm for the statistical mechan- 
ics of non-equilibrium processes. The beginning of the 
study of this phenomenon in the context of charge den- 
sity wave (CDW) dynamics was marked by the develop- 
ment of several pioneering and elegant concepts (see 0] 
for a review), in particular the description of the depin- 
ning transition in terms of critical phenomena |^. Yet 
it was hard to foresee the subsequent growth of what 
seemed to be a simple although a subtle subject into a 
fascinating multidisciplinary branch of statistical physics. 
The resurgence of interest was related to the discovery of 
high temperature superconductors (HTS) where the mo- 
tivation was driven also by the technological quest for the 
description of transport properties of HTS. 

The understanding of the remarkable effects displayed 
by driven vortex lattice involved a diversity of concepts 
drawn from various branches of contemporary physics 
ranging from polymer physics and spin glasses to non- 
linear stochastic equations and turbulence, as well as the 
invention of new concepts of non-equilibrium physics of 
disordered media (see ) . In recent years much theoret- 
ical effort has been expended to advance our knowledge 
of the driven dynamics of disordered media. Yet in spite 
of impressive achievements there remains a vast amount 
of fundamental open questions with the depth and sub- 
tleties still to be revealed. In this work we develop a reg- 
ular approach to the description of periodic media driven 
through a quenched random environment that will hope- 
fully enable to put subsequent research endeavors on a 
firm standard basis. 

Statics of disordered elastic periodic systems 

The subtle dynamic properties of dirty media are gov- 
erned by the interplay among thermal fluctuations, driv- 
ing force, and quenched disorder. To gain better insight 
into the dynamics we first discuss briefly the statics of 
weakly disordered elastic periodic systems. These include 



CDW, vortex lattices (VL), vortex arrays in Josephson 
junctions, domain walls, dislocations in solids, Wigner 
crystals, and many others. The common feature of the 
above systems is that, although the weakness of pin- 
ning suggests the purely elastic Hamiltonian as a start- 
ing point, the disorder-distorted system possesses a huge 
number of meta-stable states and the ground state is in- 
finitely degenerate. This dooms to eventual failure direct 
attacks on the asymptotic large-scale behavior based on a 
straightforward perturbation theory with respect to dis- 
order. 

The first and decisive step for an approach by statisti- 
cal physics to such systems was made in the remarkable 
work of Larkin [Q . The pioneering ideas of this work were 
later cast into the collective pinning theory 1^-0] and basi- 
cally determined the further development of the field. It 
was recognized in Q that pinning can be treated pertur- 
batively in the domain of the distorted lattice belonging 
to a single meta-stable state generated by disorder. Such 
a coherently pinned domain is called the correlated vol- 
ume and the pinning energy stored in such a domain de- 
termines the crucial characteristic of the pinned system: 
the critical depinning force. 

The key quantity characterizing the system is the de- 
gree of distortion of the elastic sy stem by disorder, the 
roughness w{r) = ([u(r) — u(0)]2), where u(r) is the dis- 
placement of a vortex from its undistorted position r. 
Within the domain w{r) < where £, is the characteris- 
tic spatial scale of variations of the random potential, the 
pinning force f has only a negligible dependence on the 
displacements u. This implies that this domain is pinned 
coherently and lies in the single valley of the effective po- 
tential landscape of the system. The roughness within 
the correlation volume grows as w^r) oc r^^ with the so- 
called wandering or roughness exponent ( that takes the 
value ^ = in the Larkin regime {d is the dimension- 
ality of the lattice). Since different Larkin domains are 
pinned independently one could conclude that pinning, 
however weak, destroys the long-range order of the lat- 
tice for d < 4. 
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As soon as the relative displacement of vortices exceeds 
the disorder correlation length, the spatial variation of 
the pinning forces becomes important. The vortices start 
to feel that they are in a random potential having many 
meta-stable states rather than in a random force field. 
Therefore the perturbative result would overestimate the 
actual roughness of the lattice. On intermediate scales 
where < w{r) < and a is the lattice constant the 
system adjusts itself to the multi-valley potential relief. 
This regime is referred to as the random manifold regime 
and the roughness exponent becomes smaller than in 
the Larkin regime. On these scales the periodicity of 
the medium is not yet of significance. The region of the 
largest scales, where lattice displacements exceed the lat- 
tice constant, < w(r), and the periodic nature comes 
into play, was first investigated by Nattermann who 
found, by the renormalization group approach, that at 
large distances pinning of the VL is equivalent to pin- 
ning of CDW and that displacements grow only logarith- 
mically, w{r) cx ln(r). The above results were confirmed 
later by variational replica approaches JTo|-p^ . 

Structures with logarithmic roughness are well known 
in the physics of surfaces and 2D crystals. The logarith- 
mic roughness implies that the system retains its peri- 
odic character and Bragg peaks in the structure factor 
S'(q), the singularities however have an algebraic charac- 
ter, S{q) oc |q — Ql^'', (Q is a reciprocal lattice vector) 
rather than the (5-function-like character as in ordinary 
crystals or Lorentzian character in liquids [p^-p^. This 
algebraic behavior is a characteristic feature of quasi- 
long-range crystalline order. 

The roughness of the lattice structure implies a rugged- 
ness of the potential landscape of the system and the ex- 
istence of infinitely high barriers separating the different 
meta-stable states, which is the characteristic feature of 
glassy systems (see This was realized in a seminal 
work by M.P.A. Fisher who identified the VL dis- 
torted by disorder as a glassy structure and called it the 
vortex glass. It is important to stress that the derivation 
of the above features was based on the elastic nature of 
disordered lattices. 

The stability of the elastic vortex glass with respect 
to the formation of topological defects (dislocations) was 
questioned by Fisher, Fisher and Huse [ISj who stated 
that dislocations are to be generated at the scales where 
the roughness becomes of order the vortex spacing and 
that therefore the elastic description of the vortex glass 
fails. In spite of the fact that the correctness of the ar- 
guments of [ p^ was questioned in its turn (the energy 
of the dislocation formation was underestimated and the 
logarithmic smoothness of the lattice on large scales was 
overlooked), the image of the vortex glass as a disloca- 
tion saturated medium became widespread. Arguments 
demonstrating the self-consistency of the elastic vortex 
glass approach (as long as the disorder is weak enough) 
were presented in ||l^,|l9|,^ . Thus the existence of a weak 
disorder-induced elastic vortex glass free of topological 
defects can be considered as well established. Recall- 



ing that a logarithmically rough medium shows algebraic 
Bragg peaks, Giamarchi and Le Doussal Q proposed 
to call the vortex lattice deformed by disorder a "Bragg 
glass." This name gained popularity among the spe- 
cialists and replaced the somewhat compromised "vor- 
tex glass" in their technical jargon. The latter terms is 
now reserved for the topologically disordered vortex solid 
phase. 

Dynamics of disordered elastic periodic systems 

The main feature of the driven dynamics in a random 
environment is the existence of the pinning threshold: at 
zero temperature the system remains pinned if the drive 
does not exceed a certain threshold value fc, the critical 
pinning force, and slides if f > fc. At finite temper- 
atures the sharp transition is rounded by thermal fluc- 
tuations and is not very well defined. Nevertheless, one 
can keep the notion of a temperature dependent critical 
force fc{T) as the force separating the pinning dominated 
regime with slow thermally activated dynamics or creep 
at low force f ^ fc from the fast sliding one at f ^ fc- 

Motion of the system in the creep regime occurs via 
thermally activated jumps over energy barriers separat- 
ing different meta-stable states. The size of the typical 
energy barrier as a function of the driving force can be 
related by scaling arguments to structural features of the 
system, in particular its roughness. This was done in l2l| ] 
for driven elastic manifolds and it was found that the bar- 
riers controlling the motion diverge algebraically at small 
driving forces as £/(/) cx 1//^, where ^= (2^— 1)/(2— ^). 
The approach of pi| was extended from continuous me- 
dia to the creep of vortex lattices in |^J^ . The divergence 
of the activation barriers implies a non-linear response 
of the system to small forces and leads to the identifi- 
cation of the low-temperature vortex state as a glassy 
phase, since such a non-linear response is a hallmark of 
the glassy system. 

The understanding of the critical behavior at the de- 
pinning threshold has seen remarkable progress p^-p^] 
since the first work by D.S. Fisher [|. 

The high-velocity sliding of the periodic systems 
was long considered as the most "easy-to-understand" 
regime. The outburst of the interest in the flow regime 
at large driving forces well above the depinning thresh- 
old was triggered by the prediction |^ of dynamic phase 
transitions between plastic sliding in the nearest vicinity 
of depinning and coherent motion of the crystalline struc- 
ture at high drives. Already early experiments p^ ] have 
shown that a moving vortex lattice has a higher degree 
of crystalline order than a pinned vortex lattice. These 
studies have been refined recently |^ 28| to identify the 
different dynamical regimes. 

The non-equilibrium phase transition predicted in Ref. 
p5[ is expected to occur in systems with sufficiently 
strong disorder, where depinning is accompanied by the 
massive production of topological defects |^ . The struc- 
tural order improves at large driving forces because the 
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system experiences disorder forces that are temporally 
fluctuating in the moving frame. In this sense the effect 
of disorder resembles the thermal noise of a heat bath. 
However, this comparison does not carry too far, since 
the quenched nature of the disorder still implies infinite- 
ranged spatio-temporal correlations of the effective force 
in the moving frame. 

Balents and Fisher |Q used scaling arguments to ex- 
tend the concept of the non-equilibrium "freezing" tran- 
sition of Ref. to CDW and have shown that true 
long-range order is restored at large velocities only in 
d > 3. Thus the critical dimension is reduced by one 
in comparison to the static case. For d < 3 the CDW 
phase, which is the analog of the vortex displacement, is 
still rough even even at the largest velocities, whereas a 
temporal order resulting in narrow-band noise may still 
persist. 

Giamarchi and Le Doussal addressed the question 
of spatial large-scale order in the driven VL. They ar- 
gued that the transverse periodicity of the system leads 
to glassy features of the driven phase. It was argued, in 
particular, that on the largest scales the driven lattice 
retains a logarithmic roughness in the directions trans- 
verse to the velocity for 2 < d < 3. In the direction 
parallel to the velocity the roughness is expected to be 
even algebraic as for driven CDW |^,^ in contrast to 
the static case, where roughness is logarithmic. One of 
the manifestations of the transverse glassiness suggested 
according to Ref. |^l| would be the existence of diverging 
transverse barriers. This brings to mind the early obser- 
vation of Schmid and Hanger who have performed 
a lowest-order perturbative calculation for the pinning 
force and noticed a discontinuity in the transverse I-V 
characteristic in a sliding state. 

Notice that in comparison to the equilibrium situa- 
tion in the absence of driving forces, the approach to 
the physics at large velocities is even more intricate be- 
cause of the non- equilibrium nature of the driven state. 
Although important predictions about this state have 
already been formulated, a systematic approach is still 
lacking. 

In this paper we develop such a systematic approach 
to the driven dynamics of dirty periodic media on the 
basis of the Martin-Siggia-Rose (MSR) formalism. The 
MSR formalism provides a powerful tool to access the 
largest scales and to treat the immediate vicinity of the 
depinning transition. Although the complete description 
should include the derivation of renormalization group 
equations, we will show that a number of important con- 
clusions concerning the properties of the driven state can 
be successfully achieved even within the framework of 
the dynamic perturbative approach. For large velocities 
V the small parameter for the expansion is Ao/(77^i;^^'^+^) 
where Aq is the d-dimensional spatial integral of the po- 
tential correlator of a width ^, and rj is the friction coef- 
ficient. 

Summary of results 



Using a coarse-graining procedure for the dynamics of 
periodic media, we find that their large-scale behavior is 
governed by an effective equation of motion 
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Renormalized parameters carry an asterisk to distinguish 
them later on from the unrenormalized (bare) ones. Un- 
der this procedure the parameters become anisotropic 
since the velocity identifies a particular spatial direction. 

All components of the friction coefhcient 77* are found 
to be finite. Therefore glassy features, which in general 
appear as a divergence of such coefficients, are absent. 
The elastic dispersion that reads -D*^(q) = —ix*apa^a + 
'^afSablalb in Fouricr space, includes besides elastic con- 
stants K also stress terms x after coarsening. 

Due to pinning and dissipative effects on spatial scales 
smaller than the coarse-grained cutoff scale (with mo- 
menta larger than A<), the velocity-dependence of the 
friction force F^'^* becomes non-linear. 

The fourth term in Eq. (|l|) is a Kardar-Parisi-Zhang 
(KPZ) non-linearity A*, which is absent in the bare 
dynamics, is generated by disorder. It is an anisotopic 
generalization of the term iA(VM)^ familiar from surface 
growth and Burger's equation. 

The pinning force /*, which was simply the gradient 
of a random potential in the bare case, acquires a more 
virulent random-force character with a correlator 
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^0 



for d > 3 , 



^{aA<f-^ ford<3 



(2a) 
(2b) 
(2c) 
(2d) 



for large velocities v ^ c£,/{rja?) with a typical elastic 
constant c. One sees that the variance of the compo- 
nents perpendicular to the velocity diverges in d < 3 as 
the coarse-graining cutoff A* — > 0. 

The effective thermal noise ^* describes in general an 
effective heat bath with a temperature that is increased 
due to shaking effects exerted by the pinning on the 
medium. 



(3a) 
(3b) 
(3c) 
(3d) 



Kp^ / (C(r,t)a(O,0)) , 



yv 



Ao 



Q^Aq 



1? for d > 3 , 
{aA<f-^-& for d < 3 
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In the non-equilibrium case this effective "temperature" 
is defined as the integral over the correlator of ^* , that 
can be distinguished from pinning forces by the temporal 
decay of its correlations. The behavior of d* and 0* is 
very similar, they show the same type of divergence for 
^0. Since this divergence comes from small mo- 
menta, it is a measure for the strong fluctuations of the 
medium on large scales only. 

The disorder generated stress couplings X: the KPZ 
non-linearities and the random forces are specific non- 
equilibrium terms that are absent in the equation of mo- 
tion before coarse-graining. 

From the effective equation of motion (|l|) the displace- 
ment fluctuations are found to roughen the medium in 
dimensions d < 3 with high anisotropy for directions par- 
allel and perpendicular to the velocity (which we choose 
to be parallel to the a;-axis). The anisotropy manifests 
itself in two distinct features of the displacement cor- 
relations: (i) The components of the displacement field 
perpendicular to v exhibit stronger fluctuations than the 
components parallel to v. (ii) The relative displacement 
increases faster in perpendicular directions than in the 
parallel direction. 

The former feature, obtained below within perturba- 
tion theory for the elastic medium, can be understood 
qualitatively already from a single-particle picture in 
analogy to the consideration that led to the notion of 
the shaking temperature, describing the disorder-induced 
increase of the effective system temperature ]25| ]. To 
this end we consider a particle moving in a disorder 
potential V with Gaussian correlations V{IV)V{0) ^ 
Ao^-'^e--^'/(2«')_ A particle starting at IL{t = 0) = 
moves with an average velocity v « F = Fx following 
an over-damped equation of motion ?7R = F — VT^(R) « 
T^v— Vl^(vt). The components of its displacement u{t) — 
R(i) —vt parallel and perpendicular to v have a variance 
that grows differently as a function of time. When the 
effect of the pinning forces is integrated over time (i.e. 
along the direction of motion of the particle), the force 
component parallel to the direction of motion is "recog- 
nized" as the gradient of a random potential, whereas the 
perpendicular components can not be distinguished from 
a true random force, since the particle does not explore 
these directions. Therefore u^^t) « ^2^2°^^ saturates for 

large times, whereas u^(t) ~ i^^'d+i \t\ grows without 
bounds, like under the influence of thermal noise. This 
implies that the shaking temperature Tgh ~ 4r ^® asso- 
ciated with the perpendicular displacement components. 

The second feature of anisotropy, a more rapid growth 
of the relative displacements w{r) in the direction per- 
pendicular to the motion (i.e. for r _L v), is related to the 
size and shape of the dynamic Larkin domain, in which 
the pinning forces act coherently on the elastic medium. 
Since has much weaker fluctuations than Uy the cor- 
relation lengths are determined by the fluctuations of Uy 
alone and are found to be 



Vc^al ^^^^ j , Xc^r]vyjc. (4) 

For weak disorder they are finite only in dimensions d < 3 
(tJc still depends logarithmically on v in d — 3) and in- 
crease for large velocities much faster parallel than per- 
pendicular to the velocity (see Fig. 

The next important question of the stability of the lat- 
tice with respect to plastic relative displacements of vor- 
tices can be captured by a phenomenological Lindemann 
criterion, that examines the fluctuations in the relative 
distance of neighboring vortices (bond length). Vortices 
neighboring in a direction parallel to v have much weaker 
fluctuations in their relative position than neighbors in 
perpendicular directions. When the relative fluctuations 
of certain bonds exceed a certain fraction of the vortex 
spacing, these bonds are expected to be broken by topo- 
logical defects. We find that the bonds in directions per- 
pendicular to the velocity have the strongest fiuctuations. 
This consideration therefore supports the suggestion that 
above a certain critical value of the velocity the VL moves 
coherently like a solid and the topological order of the 
lattice may be preserved despite of the roughness of the 
lattice in d < 3. Below this critical velocity the motion 
is plastic and vortices may move in decoupled channels. 
It is essentially the anisotropy of the Larkin domain that 
provides decoupling of flowing vortex channels. 

Fig. 1^ summarizes our view of the dynamic phase dia- 
gram for the case of sufficiently strong disorder. Starting 
from highest drives we expect a coherent motion of the 
topologically ordered phase. Upon decreasing the driv- 
ing force the fluctuations of the bonds between the neigh- 
boring vortices cause a massive production of topologi- 
cal defects at the transition from coherent to incoherent 
motion marked by the solid line. This line corresponds 
to the freezing transition of ||2^. The question concern- 
ing the nature of the plastically moving phase still re- 
mains. Our analysis suggests that there is a tendency to 
channel formation, but at this point we cannot conclude 
whether these channels remain stable upon a further de- 
crease of the applied force, and therefore dynamic melt- 
ing describes transition from the moving quasi-crystal to 
moving smectic or directly into the fluid like phase. 
The possible transition between the smectic and fluid-like 
phases is denoted by the dotted line. The lower strip be- 
low the critical current (the dashed line) corresponds to 
the pinned state where the system moves via thermally 
activated jumps between meta-stable states. 

For weak disorder a dynamic transition from the co- 
herently pinned phase to a coherently moving phase is 
possible without passing through a plastic regime. Plas- 
ticity occurs only at sufficiently high temperature and for 
small enough velocities. Since the anisotropy of the sys- 
tem decreases with decreasing velocity, the width of the 
smectic regime shrinks in that direction. 

In fact it also remains a fundamental open question to 
what extent the creep regime can be considered as co- 
herent in the sense that the topological order persists up 
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to the largest length scales. The successful description of 
this regime by collective pinning requires only the typi- 
cal distance between free topological defects to increase 
faster than the size of the largest effective barriers for 
decreasing creep velocity. In principle it is possible that 
only at strictly zero velocity the coherence of the lattice 
is restored. 

The paper is organized as follows. In Sec. p we spec- 
ify the model under consideration. In Sec. Ill the gen- 



eral perturbative approach is established in a dynamical 
formalism. A scheme for the systematic extraction of 
coarse-grained parameters that describe the physics at 
large scales is presented. These parameters are evalu- 
ated in Sec. IV and lead to our conclusions in Sec 0. 
The complexity of the problem requires a compactified 
notation that is summarized in Appendix ^ Interme- 
diate steps of our calculations are sketched in Appendix 



with the Bardeen-Stephen friction coefRcient r], the elas- 
tic force Du ~ cd^u to be specified below, a driving 
force F, a thermal noise ^(r, i), and a pinning force 
f(R) ~ —dVCR). Both the thermal and pinning forces 
are supposed to have a Gaussian distribution with zero 
average and correlations 



y(Ri)y(R2) = A(Ri2) . 



(7a) 
(7b) 
(7c) 



To make our formulas comprehensive and transparent, we 
introduce a shortened notation, where e.g. ri2 :— ri — r2 
and ti2 :— ti — 12 (see also Appendix ^ for definiteness). 
Greek indices represent components in the d± directions 
of r I . 



II. MODEL FOR DRIVEN VORTEX LATTICES 



A. Action formulation 



To be specific we introduce a model of a d-dimensional 
vortex lattice. The most common realizations are vortex 
lines in a three-dimensional superconductor (d\\ — 1 and 
d_L = 2), point vortices in a film {d\\ = and d± — 2), 
or vortex lines confined to a plane (du — 1 and d±^ — 1). 
We use a unifying description by considering vortices as 
d|| -dimensional objects that can be displaced in d±^ direc- 
tions within a d = c?|| -|- dx-dimensional superconductor. 

To every individual vortex we assign a fixed label rj^ 
that coincides with its position in a perfectly ordered 
lattice. The coordinates along the magnetic field are de- 
noted by r||. The actual position of vortex at time t 
is denoted by R(r, i), where r — (r_L,r||). For a three- 
dimensional vortex line lattice r_L is a vector in the (x, y)- 
plane and rn = Rn represents the z-coordinate, see Fig. 

i . . 

We consider a sliding state where vortices move with 

the average velocity v. Then r is viewed as the undis- 
torted vortex position in a comoving frame, whereas R is 
the actual position in a laboratory frame. To parameter- 
ize the fluctuations of the vortex lattice, we define vortex 
displacements as 



u(r,i) = R(r,i) - r - 



(5) 



The proper choice of the perfect lattice positions and of v 
guarantees that these displacements always vanish upon 
averaging over thermal fluctuations and disorder. 

We restrict ourselves to the elastic lattice, where the 
topology of the vortices is fixed and their interactions can 
be treated in the harmonic approximation. The dynam- 
ics of the vortex lattice is governed by the over-damped 
Langevin equation 



77u(r, t) = Du(r, <) + F - r^v + ^(r, t) 
+f(r + vi + u(r,t)) 



The main difficulty in solving Eq. (g) is the highly 
non-linear dependence of the pinning force on the dis- 
placements. We will treat this non-linearity by a pertur- 
bative expansion m. l/v. A convenient way to explore 
dynamics is the standard field-theoretical representation 
of Martin, Siggia, and Rose | p^ , p7t . In this formalism 
the partition function for the out-of-cquilibrium system 
is defined as 



Z := 



j Vts\u,u\ 



(A) 



V^[u,u] := [|d'^^u(q,i)d''^u(q,t) , (8) 

where the path integral is restricted to modes g < A with 
the cutoff A This scale can be related to the coher- 
ence length as "diameter" of the vortices. The auxiliary 
response field u is introduced in addition to the displace- 
ment field u. 

To every possible configuration of the fields (includ- 
ing their time-dependence) a statistical weight is 
assigned with an action A = A[u^u]. The sum over all 
weights is normalized to unity and is independent of the 
random pinning potential. Therefore a disorder averag- 
ing can be performed straightforwardly, which produces 
a translation-invariant effective field theory. We decom- 
pose the resulting action into the "pure" and the "pin- 
ning" part 

A = ylP"'<= + y^P'" , 
^P"'" = ^ l^ui • ui + mi • [yyui - (Du)i - F + r;v]| , 



(6) 



^P'" = - / Ui • *(Ri2) ■ ii2 

^ Jl2 



(9) 
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We have introduced further abbreviations in Eq. (||) 
(see Appendix]^): an integer index i stands for {Ti,ti), 
Mi = \i{vi,ti), Ri2 = ri2 + vti2 + Ui2, and the scalar 
product includes all d space components. represents 
a short-hand notation for an integration over ti, an inte- 
gration over and a summation over the vortex labels 
T±i. The latter summation includes the factor a"*^ repre- 
senting the volume per one vortex (for the usual vortex 
lattice with two displacement components = *&o/-B 
with the flux quantum $o and the magnetic induction 
B) such that the sum would become an integral in the 
continuum limit a — > 0. It is important to stress that 
although the discreteness of the lattice does not manifest 
itself in the notation, the discrete nature of the system 
is completely accounted for in our description. The pin- 
ning part contains the pinning force correlator which 
is related to the pinning potential correlator A by 



*(R) = $„0(R) := /„(R)/;3(0) = -a„9^A(R) . (10) 

By a replacement $(Ri2) ^ $(ri2, R12) the present the- 
ory can be generalized straightforwardly to other prob- 
lems like that of interfaces dynamics or random XY- 
models 

The discreteness of the vortex lattice enforces a peri- 
odicity in Fourier space 

u(q + Q,c^) = u(q,c^) := / e^'^'-^'^ u(r, t) (11) 

Jtr 

with reciprocal lattice vectors (RLV) Q. 

B. Pure part 

The elastic interaction of vortices in the harmonic ap- 
proximation is most conveniently represented in Fourier 
space, where the pure part of the action acquires the form 



iut - f -u + mt - r -u 



(12) 



The symbol u = u(q, cj), = u(— q, — cj) (the dagger 
stands for transposition of components and complex con- 
jugation of Fourier-transformed quantities), and the q- 
integration runs only over the first Brillouin-zone of the 

lattice. Note that one can always choose f ^ = f . In Eq. 
( |l^ ) we have dropped the terms linear in the response 
field. Since the average velocity is defined by the condi- 
tion that the average displacement has to vanish, these 
terms actually have to cancel each other in the absence 
of disorder. This implies F = 77V. 

It is convenient to write the propagators of the pure 
action in the normal mode representation. In the dj^ — 2- 
dimensional case these modes are longitudinal (L) or 
transverse (T) with respect to qj_. Using the projectors 



L.ap 



(q) 



gag/3 

gi 



PT,ap{ai) ■■= Sa/3 - PL,Q/3(q) (13) 



we have {p = L,T) 
f(q,c^) 



^fj,(q,c^)Pp(q) 



p 



and a similar expression for F, where 

fp(q,c^) , 

rp(q, uj) := -irjuj + Dp{q) . 

The elastic dispersion relations read 

DL{q) = ciiq\ + c^qj , 

Dt{(1) = cgegi + C44q'\\ , 



(14) 



(15a) 
(15b) 

(16a) 
(16b) 



where the elastic constants for compression cn, shear 
C66 , and tilt C44 of the vortex lattice can have additional 
implicit dependences on q pi. 

Within the partition sum ^) response- and correlation- 
functions are defined as 

G12 = Gai3{ri2,ti2) := {uairi,ti)iup{r2,t2)) , (17a) 

C12 = Cal3{ri2,tl2) := (Ua(ri,tl)u/3(r2,<2)) ■ (17b) 

These averages are to be performed with the action (||) 
that has already been averaged over disorder. There- 
fore no further disorder-averaging is required in (|l^) and 
these two-point quantities depend only on time and space 
differences. Ga/^ (r 12, ^12) = S{ua{ri,ti)) /d^f3{r2,t2) de- 
scribes the response of the displacement u to the thermal 
noise ^ or some additional external force acting on the 
vortices and is therefore causal, i.e. G(r,t) — for t < 

ii- 

In the sequel the average squared displacement of vor- 
tices 



Co,Q/3 := (ua(r, i)u0(r,t)) 



(18) 



will play a particular role, as well as the difference corre- 
lation 

Wal3{ri2,ti2) := ^{[Uia - W2a][wi/3 - U2p]) 



— Coaf 



-jCafi{vi2,ti2) 



^C^f3ir2ut2l) . (19) 



For the general pure action ([l^ ) the response and cor- 
relation functions are related to the propagators by 

Giq,u)^T-\q,u) , (20a) 
C(q, uj) ^ G(q, oj) ■ f(q, uj) ■ G'' {q, uj) . (20b) 

The last equation implies = C in agreement with the 
definition (|l|b). 

In the special case ([l5|) we have explicitly 



Gp(q,cj) 
Cp{q,uj) 



1 



1 



rp(q,w) 
fp(q,t^) 



-irjuj ■ 



(21a) 
(21b) 



The decomposition (| 
these matrices. 



rp(q,c.)|2 ,^2^2 + ^2 (q)- 

into normal modes holds also for 
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III. PERTURBATION THEORY: GENERAL 
SCHEME 

In this Section we perform a perturbative coarse- 
graining for the dynamics of the vortex lattice. We find 
an eS'ective dynamics for the fields with wave-vectors 
q < A"^ below a reduced cutoff A< < A by averaging 
over all modes with wave-vectors A^ < (? < A. This av- 
eraging is performed by integrating out these modes in 
the partition sum (H). Due to the high non-linearity of 
the pinning part of the action this integration cannot be 
performed exactly. We restrict our perturbative analysis 
to the first and second order in the pinning action. 

As usual we separate the modes with wave vectors be- 
low and above the new cutoff, 

u = u< + u> , (22a) 
u<{r,t)= f e^'>'-u(q,i) , (22b) 

Jq<A< 

u>(r,i)- / e^'''-u(q,<) . (22c) 

JA<<q<A 

The response field u is treated analogously. The integra- 
tion over the modes u> and leads then to an effective 
action for the modes "<". The effective action can be 
represented via a cumulant expansion which is 

A''^ = A + A^^^ + ^(2) + . . . , (23a) 
^(i)[u<,u<] = (^p'"[m< +{t>,u< +u>])> 

-AP"'[u<,u<] , (23b) 

AP"'[u< +u>,u< +u>])> (23c) 

to the second order perturbation theory in ^p'". The av- 
eraging is performed over the modes and weighted 
with the pure action ^p"''''[m>, We denote the sec- 
ond cumulant by {A,A)c ■= {A'^) - {A^. 

Now we turn to a detailed analysis of the above cor- 
rections and derive new couplings in the effective action 
for the large-scale modes. 

A. First order 

Using Fourier-transforms of the disorder correlator 

^Pi" = i /" Qi . *(k) ■ U2 eJ^ '^'' (24) 

2 Jl2k 

we shift all the dependences on the displacements to the 
exponential. Applying Wick's theorem to the first order 
correction given in Eq. (|2^) we get 

A'^'\u<,u<] = / g.k.R<-k.W>.k 
Jl2k 

X |ifi< . *(k) . u< + u< • *(k) . G> t . k| . (25) 



The terms ~ GfjGj"]^ vanish due to causality. In addition 
terms ~ G-^^ etc. vanish in the Ito calculus. As a result 
all terms that appear in the correction to the action con- 
tain at least one response field. The superscripts in 
Eq. ( p5| ) stands to remind that = ri2+vti2-|-u^2 ^-^id 
that and arises from the averaging over modes 
">". To keep the notation simple, we will hereafter drop 
these superscripts. 

The obtained correction to the action has a more com- 
plicated functional dependence on the fields than the 
original action A. From we extract not only correc- 
tions to the parameters of the original A, but also new 
couplings like a Kardar-Parisi-Zhang (KPZ) nonlinearity 



1. Friction force 



The conventional approach to the description of driven 
disordered dynamics is to fix the external drive F and 
evaluate the resulting drift velocity v = v(F) as a re- 
sponse thereof. However, since in the present formulation 
F appears only at one position in the action, whereas v 
appears at many positions, we find it more convenient 
technically to treat the average velocity as given and to 
derive the force F = F(v) necessary to maintain this ve- 
locity. The fact that v is indeed the average vortex veloc- 
ity is expressed by the condition (u(q — 0,u} ~ 0)) = 0. 
In the presence of disorder this requires that the effective 
action for large-scale modes has no terms linear in u, i.e. 
no terms of order 0{u^u'^) in the fields. In the absence 
of disorder this immediately yields that the driving force 
is compensated by the friction force F^'' — r/v. 

The first-order correction to the friction force is ex- 
tracted from order 0{u^u'^) of (p5|): 

A'ilo = / lii • *(k) • 0^21 • k e*-(r-+v*i2)-k-w,..k _ 

Jl2k 

(26) 

Comparing this contribution to the term of order 
0{u^u'^) in the original action we identify the first or- 
der perturbative correction to the pinning force 

_pfr(l) ^ _j /" g4k-(r+vt)-k-W(r,t)-k 
Jtrk 

x$„^(k)G^^(r,i)fc^ . (27) 

Consequently the transport characteristics of the super- 
conductor is given by F^'^(v) = rjv + F^'^'^^'> {-v). Since F^^^ 
is in general a non- linear function of v, the characteristics 
is no longer linear over the whole current range. 
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2. Friction coefficient, elastic dispersion 

The friction coefficient 77 and the elastic dispersion re- 
lations D parameterize the propagator F in the original 
pure action ( p^ ) in the order 0{v}u^). Therefore, correc- 
tions to these parameters are extracted from (|2^) in the 
same order, 



„jk-(ri2+vti2)-k-Wi2-k 



12k 



X [{ii . *(k) . Gt2i ■ k] [zk . (ui - U2)] . (28) 



.(1) 



1 



^ik-(r+vt)-k-W(r,t)-k 



trk 



(35) 

Comparing Eqs. ( p7| ) and (|3^) one finds straightfor- 
wardly the following useful relation: 



(36) 



In the first order perturbation theory the friction coeffi- 
cient coincides with the differential resistivity. 



The integral kernel now has a finite width in (ri2,ti2)- 
However, this width is smaller than that of the response 
function, which decays on a characteristic scale in 
space and on a scale rj/D{K^) in time. 

Aiming at the physics at scales much larger than the 
width of the kernel we consider the displacement fields 
as nearly constant and approximate 



Ui - U2 W tl29tUi -I- ri2adaUi 
--^ri2ari2bdadbUi - 



(29) 



(Latin indices also include directions parallel to the vor- 
tex lines.) Inserting ( p9| ) in ( p8| ) we find 



(30) 



where the correction to the response propagator can be 
written as 



r 



(1) - -i„^^^cj + D^^^(a) 

a[3 — ^^aff^ + ^a/3\^) 



with an elastic dispersion 



(31) 



(32) 



Therein new coefficients x appear. They describe forces 
arising from a direct coupling to lattice stresses. In ad- 
dition, the coarse-graining generates elastic constants k 
with reduced symmetries compared to the original con- 
stants. 

The correction to the friction coefficient is explicitly 



'trk 



ik-(r+vt)-k-W(r,t)-k 

xt$Q^(k)G5^(r,t)fc5fc/3 . 



(33) 



In general this tensor is non-diagonal and gives rise to 
Hall effects. The stress coefficients are 



ik-(r+vt)-k-W(r,t)-k 



trk 



xra<S>a^i'k)G5^ir,t)kskp (34) 
and the elastic couplings are corrected by 



3. The KPZ term 
Action (ESl) contains a further contribution in order 

/((I) [ „Jk-(ri2+vti2)-k-Wi2-k 

^ Jl2k 

X [ui • *(k) • 6^21 • k] [k • (ui - U2)]2 . (37) 

As before, we use the approximation (|29| ) which leads 
to (we omit other terms that are less relevant on large 
scales) 



(1) 

uuu 



(38) 



with 



A^'^ --II e 

/trk 



ik-(r+vt)-k-W(r,t)-k 



xrarfcfc^fc^fce$Q5(k)Ge5(r,t) . (39) 

This means that in the driven state disorder induces a 
KPZ term. 



4- Disorder correlator 
So far we have not considered 0{u'^) of (p5|): 
4'^ = 1 [ ["1 • *(k) ■ U2]e^'"^^^ [e-k Wi2 k _ _ 

^ Jl2k 

(40) 

Comparing the functional form of this action to ( |2^ ) 
we identify a correction ^'•^''(k) to the disorder corre- 
lator from the persistent part of the kernel, i.e. the 
part that is present also for |ti2l — » oo. Using 
lim|tj2|->oo W(ri2,ti2) = Co wc find 

*(i)(k) = ^{k)[e-^-^»-^ - 1] . (41) 

Remarkably, is independent of velocity and vanishes 
in a perfectly ordered lattice with Cq = 0. In the general 
case with finite Co the correction represents a smearing- 
out of disorder by the vortex fluctuations. 
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5. Temperature 



1. Random force 



The remaining non-persistent part of (|40|), which is not 
taken into account by the disorder correction, is 



^ ■/12k 



ikRi2 



g-k-Wi2-k _ g-k-Co-k-| 



X[ui •*(]<)• U2] 



(42) 



Now this integrand vanishes for |ti2| — > 00 and is also lo- 
cal in ri2. Assuming again that the width of this kernel is 
small compared to the scales of variation of u and u, one 
may neglect Ui — U2 (as zeroth-order of approximation 
(v%) and approximate 



^V-'^^i^ u(r,i).i?«.u(r,t) 



Herein 



Itrk 



k-(r+vi) 



,-k-W(r,t)-k 



-k Cn k 



(43) 



*a/3(k) 



(44) 



is the correlator of the effective thermal noise (non- 
persistent shaking forces). Note that -d^^^ = for a per- 
fectly ordered lattice with Cq = 0. 



B. Second order 

The second-order correction A^'^^ to the action has to 
be calculated according to Eq. (p3|c). The result con- 
tains a large number of terms, and the full expression 
is not displayed here. One can easily see, however, that 
A'--^^ contains terms of orders 0{u^), 0{u'^), and 0{u^). 
From 0{u^) corrections to the pinning force and to T 
can be extracted. From 0{u^) one derives corrections to 
the disorder correlator and to T. Eventually, new types 
of couplings appear in 0{u^) that represent higher-order 
cumulants of the disorder, i.e. deviations from a Gaus- 
sian distribution. 

The subsequent analysis is restricted to the evaluation 
of the correction to the disorder correlator, in order to 
demonstrate the presence of effective random forces in 
the coarse-grained equation of motion. The second order 
correction to the disorder correlator is again ob- 

tained by identifying the contributions to 0{u^) in A^^"^ 
that represent persistent (quenched) forces in the labora- 
tory frame in contrast to temporarily fluctuating forces 
that contribute to the effective thermal noise. 

Due to the complexity of the involved expressions the 
technical details of this procedure are deferred to Ap- 
pendix H. 



Now we determine the coarse-grained disorder correla- 
tor. The second-order correction is extracted from A^^'^ 
in 0{u^). We identify ^'■^-'(k) as the kernel where slowly 
varying displacements and response fields enter in exactly 
the same combination as they appeared in the original 
pinning action (p^. Therefore $^^-'(k) represents the 
correlator of forces that are stationary in the laboratory 
frame. The force experienced by a vortex moving in the 
laboratory frame is nevertheless fluctuating in time. 

The calculation of this correlator is performed in Ap- 
pendix IB^ and leads to the somewhat involved expression 
for *(2^) given in Eq. (||) . From the large scale behav- 
ior (k 0) 



$a/3(k) = -I- (j)al3fk^ + -^(j^aPjSkjks + . . . (45) 



one identifies the random force correlator (t>a/3- This con- 



tribution as well as the second term in (45) emerge only 



in the driven state and in the presence of disorder. The 
bare random potential contributes only to the coefficient 

</>Q/375 = 2Saj6f3sA{k = 0). 

Eq. (0) gives explicitly 



-(2) 



,-k-[W(ri,ti)+W(r2,i2)]-k+jk-R<'' 



12k 



X {a>„0(k)$^5(k)e^'^-^="' - $„^(k)<i>;35(k)e-^''-^^'" } (46) 
with rI"^ := r, +wti. 



2. Other terms 



One can also extract coefficients (fiap-y and (jjap-yS of the 
force correlator from (||) . However on large length scales 
the corresponding terms in the action are less relevant 
than the random force, and therefore we do not present 
them here. 

A number of other terms appear in the second order of 
perturbation theory introducing in particular new types 
of disorder. For example, taking into account the gradi- 
ent term in the expansion (^ ) one finds a random cor- 
rection to the amplitude of the KPZ nonlinearity as sug- 
gested by Krug . Also the second-order corrections to 
F and to the propagators appear. Again, these correc- 
tions have a complicated form and are not given here. For 
weak disorder the second-order corrections are expected 
to be small compared to the first-order corrections. 
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C. Fluctuation-Dissipation Tlieorem 

In the absence of disorder the system obeys the 
Fluctuation-Dissipation Theorem (FDT) ^] 

- 2i77wf„/3(q, w) = i9[r„/3(q, i^) - r/3a(-q, -lo)] , (47) 



which imphes (for i 7^ 0) 



(48) 



The FDT holds even in the driven state due to the 
Galilei-invariance (the response- and correlation func- 
tions under consideration are defined in the comoving 
frame) . 

In the presence of disorder the validity of the FDT 
would require 



- 2z77c.ri^j(q,c.) = t?[rl^^(q,c.) - r^^^(-q, -c.)] . (49) 

At zero velocity this relation is satisfied even in the 
presence of disorder, because the terms x vanish and 
Tyi?!^^^ — 'di^'^p . However, in the driven dirty system the 
FDT is violated. The most obvious reason is the pres- 
ence of the stress couplings generated by disorder in the 
driven system, also one sees immediately that 'd'^^^ and 
rj^^^ no longer satisfy the aforementioned relation. 



(1) 



IV. EVALUATION OF PERTURBATION 
THEORY 

In the previous section we have derived the effective 
equation of motion (|l|) for the coarse-grained displace- 
ment. The renormalized parameters (with superscript 
*) are composed of the original values plus perturbative 
corrections, 



Dip 



Dal3 - 



, „(1) 
'lap ' 



A! 



(C(ri,ti)?;(r2,t2)) 

Kp 

/*(Ri)/S(R2) 



= A 



(1) 



(1) 



(1) 

Q/3 



(2) 

Q/3 



Now we evaluate the various couplings that appeared 
upon the coarse-graining procedure. We start wit h the 
first-order perturbative corrections in subsections IV A 



IV G and then derive the random-force term as the 



second-order correction in subsection tV H 



To simplify the evaluation of the general perturbative 
results we consider disorder with an isotropic correlation 
length ^ 



A(k) = Aoe- 



(50) 



where we assume ^ ^ a, which is typical for supercon- 
ductors. In addition, we assume the elastic constants to 
be uniform, i.e. Dap{ci) = Sapcq^ whenever explicit ex- 
pressions involving elastic constants are given. This will 
not change our results qualitatively. 

As shown by Schmid and Hauger the lattice prefers 
to move along the principal symmetry axes. Hereafter we 
restrict the analysis to the situation where the velocity 
is parallel to a high-symmetry direction of the lattice, 
which we choose to be the x-axis. 



A. Random potential 

We start with the discussion of the coarse-grained 
disorder correlator. The first-order correction (|4l| ) 
preserves the random-potential nature of the original 
disorder. Coarse-graining smears out the correlation 
length of the disorder over the typical vortex displace- 



ment since $a/3(k) + $i^j(k) = e-^-^°'^<^apO^) = 

kakpAoe~^''-^'''^~^^ ^ . This means that the disorder 
correlation lengths is described by the matrix 



Sab £,'^Sab + Coab 



(51) 



The correction Goab diverges in dimensions d < 2 for fi- 
nite temperatures, if the lower cutoff A^ is sent to zero. 
In this case it is therefore possible that weak disorder is 
irrelevant for the large-scale properties of the vortex lat- 
tice. As long as A< > finite this correction has only a 
quantitative effect in all dimensions. 

Focusing on c? > 2 we will consider in what follows 
mainly the case of zero temperature. In most perturba- 
tive expressions the exponential factors involving matri- 
ces Co or W can be ignored since they always enter in a 
combination with the disorder correlator and modify the 
correlation lengths only quantitatively. The only excep- 
tion is provided by the correction to the temperature. 



B. Temperature 

We evaluate the effective temperature from Eq. (|44|), 
where we had observed already that the correction van- 
ishes for the perfectly ordered lattice at T = 0, where 
W = Co = 0. At finite temperatures one expects a 
positive correction, since in general W < Co and the 
difference between the exponentials in Eq. (^) will no 
longer vanish. Considering low temperatures, one may 
linearize the exponentials and obtains 



(1) 

a/3 



(52) 
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where we implicitly decompose k = Q + q into a recip- 
rocal lattice vector (RLV) Q and a vector q within the 
first Brillouin zone. With a random-potential disorder 
(^) one immediately finds in the limit of large velocities 



^(1) 

XX 



^(1) 
yy 



^(1) 
yy 



Ao 



Q^Aq 



i9 for d > 3 , 
(aA<)''-3^ for d < 3 . 



(53a) 
(53b) 
(53c) 



We have dropped numerical factors of order unity. For 
large velocities the result for -d^xx is independent of the 
elastic constants and finite for A< — * 0. 

However the integrand of i}[^J has poles near k w Q 
giving rise to the divergence in d < 3 (the divergence 
is logarithmic in d = 3) when the coarse-graining cutoff 
A< is sent to zero. This divergence is the first indica- 
tion that perturbation theory is insufficient for the true 
large-scale description of the drifting vortex lattice even 
for arbitrarily large drive. 

Since the correlator '&y^J decays even for finite A< with 

a smaller power in v than "d^x , the effective noise is 
stronger in the directions perpendicular to the velocity. 
Since only i^l^J exhibits this divergence, there is a funda- 
mental difference between vortices with dj^ > 1 displace- 
ment components and CDW with d^ = 1. 

This divergence therefore does not imply that fluctu- 
ations, which have on small scales an amplitude propor- 
tional to temperature in the pure case, would now di- 
verge. Since this divergence of the effective temperature 
arises from the large-scale response of the elastic medium 
to the pinning force, it suggests only that on large scales 
the medium is much more rough than in the pure case. 
However, the roughness does not depend on the effec- 
tive temperature only, but also on the effective friction 
coefficients and elastic constants, which we thus have to 
evaluate before we can address the roughness. 



C. Friction force 



where purely numerical prefactors of order unity have 
been dropped. This expression is independent of the elas- 
tic constant for large velocities. In this regime vortices 
respond dynamically like individual particles. The typi- 
cal force scale Fq = Aq / ^'^+2 is set by the spatial average 
of the pinning force. 

Eq. ( |55| ) applies only to velocities parallel to a prin- 
cipal lattice axis. For other directions the force and ve- 
locity will no longer be parallel to each other [Q: the 
velocity will deviate from the force in the direction of 
the closest principal lattice axis. Thus disorder induces a 
Hall-effect. Instead of evaluating the friction force for ar- 
bitrary directions, we will turn to the friction coefficients 
that describe the same effect in a differential form. There 
the presence of the Hall-effect shows up as anisotropy of 
the coefficient matrix. 

At small velocities Eq. (|5^ ) gives Ohmic behavior, 

Fx^^"^ ~ V, only for d > 4. In this case the poles of the 
integrand at k = Q with Q • v = are integrable. This is 
no longer true in d < 4 and the perturbation theory gives 
sub-Ohmic transport Fx^^'^ ~ yi^-"^)/^^ in particular in 
d — 2 the effective friction force is finite for small veloci- 
ties |3^ . The sub-Ohmic behavior of the effective friction 
force reflects the fact that the vortex lattice forms a glass 
below four dimensions and at zero velocity M . 



D. Friction coefficients 

In order to examine to what extent glassy features per- 
sist in the driven state at finite velocities, it is necessary 
to examine the friction coefficients rjafs- Since these co- 
efficients describe the dynamical response of the driven 
vortex lattice, the glassy features, which are in general 
associated with divergent relaxation times, must mani- 
fest themselves as divergences in the friction coefficient. 

Combining (54) with (|3q) we find 



77 / fc„fc^fciA(k) 



rj'^v^k^ 



[rf-v'^kx 



(56) 



Note that only the diagonal components of the friction 
coefficient do not vanish due to reflection symmetries. In 
the limit of large velocities parallel to a principal lattice 
axis (x-axis) one finds: 



At T = and uniform elasticity the evaluation of Eq. 
(|^) yields 



(54) 



The relative correction to the friction force closely related 
to the relative correction (^3|) of the effective tempera- 
ture. We find explicitly for v » ^c/yya^ 



Ar 



(55) 



„(1) 

'I XX 



(1) 
lyy 



Ao 



aAo 



77 



(57a) 
(57b) 



The main contributions to rjy]^ come from the vicinity 
k « Q of RLV with Q • v = 0, whereas for rji 
Q • V ^ contribute. We find that \r]y]^\ > 
we have evaluated Eq. (|5^) in the limit ^ ^ a. Oth- 
erwise both corrections would be of the same order of 
magnitude. 



■ji]^ all terms 
Iryi^l since 
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Since ( |55| ) has already been specified to velocities along 
a principal lattice axis, only ryi^ can be derived from 
there using the differential relation (|36|). Both friction 
coefficients are again independent of elastic constants. 

In their pioneering work Schmid and Hanger have 
discussed a discontinuity in the relation between the 
transverse force and velocity [see their Eq. (27) and dis- 
cussion thereafter]. However, as they state, this discon- 
tinuity is an artifact due to a neglect of vectors q. Such 
a discontinuity, if real, should have appeared as a diver- 
gence in liyy, which is actually absent [Q. 

Glassy features in the dynamical response normal to 
the velocity, suggested by Giamarchi and Le Doussal , 
should manifest themselves as a divergence of the friction 
coefficient r^yy . At large drive this divergence can emerge 
only in higher orders of perturbation theory. 

At zero velocity a divergence rj^^p ~ A^'' * is appar- 
ent for d < 4. It arises from the poles of the integrand in 

(H) at k = Q and implies yy^^j - v^'^-'^^/^Saf} for A< = 0. 
This again signals the glassiness at v =0 for d < 4. 



E. Stress coefficients 



For T = and uniform elasticity Eq. (| 



c^qaq 



reduces to 
(58) 



We again restrict our consideration to the case where the 
velocity is parallel to a principal lattice axis. For several 
combinations of indices the integrand is odd under an 
inversion ka —ka- Therefore we find in particular 



(1) 



A-yxx 



Axyx 



A^xxy 



■^yyy 







(59) 



(We would like to remind that Greek indices run only 
over the directions perpendicular to the vortex lines, i.e. 
X and y in the usual 3D configuration, whereas the Latin 
indices also include directions parallel to the vortex lines, 
i.e. z in the usual 3D configuration.) But for finite veloc- 
ity there are also non- vanishing components which decay 
in the limit v oo like 



A.XXX 



Aoc2 



■^xyy -^yxy 



Aoc 



Ayyx 



Aqc 



^d+3j^2y2 



(60a) 
(60b) 
(60c) 



We have assumed ^ ^ a. Then Xyyx is the largest among 
these coefficients. 

The stress coefficients have a simple physical meaning. 
They express the tendency of vortices to form a homo- 
geneously moving system. In particular the coefficients 
Xxxx > and Xyyx > imply that a vortex experiences a 



dynamical force that makes it follow the motion ("foot- 
steps") of the precedent vortex, see Fig. |[ These cou- 
plings favor the formation of channels. 

The stress coefficients are generated only in a non- 
equilibrium driven state and in the presence of disorder. 
At small velocity and for A^ ~ they vanish ~ v only in 
dimensions d > 4, where the integral in ( |5^ ) is finite. For 
d < 4 the stress coefficients vanish sub-linearly. Lowest- 
order perturbation theory gives e.g. Xxxx ~ u(''-2)/2^ 
Thus these coefficients diverge in d < 2 only. 



F. Elastic constants 



Following the same scheme a correction to the elastic 
constants 



^a(3ab 



A(k){,5a6 



rfv'^kf 



-4c qaqbq' 



i-q^v'^kl 



[ifv'^kl+c^q^Y 
} (61) 



2 4 

c^q 



can be obtained from Eq. (35) 



For symmetry reasons again all coefficients k^^^^^i, van- 
ish where indices different from x appear an odd number 
of times. Nevertheless k^^^^i, will no longer be propor- 
tional to SapSab and the correction reduces the symmetry 
of the original elastic interaction. 

The first term in ( |6l| ) represents a correction of the 
elastic constant proportional to the correction of the fric- 
tion coefficients (56). This contribution is finite for all 
dimensions and of the order v~'^ for large v. The sec- 
ond term in mU) can be shown to be finite and of order 



^ V . Therefore, at large velocities, it can be neglected 
in comparison to the first term. 

We find thus in the leading order ~ v^^ 



A, 







,(1) ^ 

^xxab ~ ^d+2j^2y2 



Sab C , 



(1) 

yyab 



Sab C 



(62a) 
(62b) 



Hence the elastic stiffness of the lattice for displacements 
parallel to the velocity is reduced, whereas the stiffness for 
displacements perpendicular to the velocity is increased. 
The latter effect can be interpreted as a tendency to form 
channels and an increased energy cost for transverse dis- 
placements. 

For small velocities the corrections to the elastic con- 



stants are found to diverge k*-^-* 
tion coefficients. 



/d-4)/2 ijj^g ^Yie fric- 
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G. KPZ nonlinearity 



As before we obtain from Eq. (P9| ) 



Acqv 



A(k) 



SabCq'^ 



+'2cqaqb 



[if'v'^k'^ + c?q'^Y 



[rfv^kl + c^q'^Y 
} ■ (63) 



Many terms vanish due to symmetry. All coefficients 
that are related by permutations among Greek or among 
Latin indices are identical. 
We find in particular 



^(l) ^ g'Ap 
^xxxab ~ ^2^d+2jj3y3 



^ab , 



X 



(1) _ cAq 



3ah , 



xyyab ^d+3j^2y2 

Ad) 

^yyyyx ^^d+s^^a^a 



(64a) 
(64b) 
(64c) 



All couplings assume finite values and decay at least 
^ v^"^ in the limit of large driving. 

For CDW, which are included in our analysis by spe- 
cializing all Greek indices to x, Chen et al. j45j have found 

'^^xxxab in contrast to our result (|64|a). Having no 

access to their derivation, we were not able to pinpoint 
the origin of the disagreement. 

In the limit of small velocities the KPZ terms coeffi- 
cients vanish ~ v only as long as the integral in ( |63|) is 
finite, i.e. for d > 6. For d < 6 the coefficients scale like 
A ^ ti('^^4)/2 Qj-^g ^^^^ fhai these coefficients diverge 
even for d < A\ 



H. Random force 

At T = and uniform elasticity the evaluation of the 
random force correlator Eq. (Efl) yields 



b^0-2JJ^kpk^A (k) 



(65) 



This expression was recently given in p3^ where it appar- 
ently was found within a RG framework. Here we obtain 
it as a result of a straightforward perturbation theory. 
One easily calculates for v ^ £,c/rja^ 



,(2) , 


A^ 


XX 








yy 




,(2) , 


. «'A2 


yy 





for (i > 3 , 



^(aA<)''-3 ford<3 



(66a) 
(66b) 
(66c) 



(2) 

For large velocities the result for (frxx is independent of 
the elastic constants and finite for A< 0. This random 



force is the analogon of the random mobility of driven 
interfaces with phase disorder [^9| . A similar force for 
CDW was previously obtained in [^,^ and predicted 
for vortex lattices also in [^,0 . 

Comparing this result to the noise correlator (^3|), we 
find exactly the same type of divergences in the trans- 
verse components in the limit A^ — > 0. 

At small velocities the random force correlator ( |65| ) 
vanishes like in d > 8 and ~ yi^--'^)/'^ in lower di- 
mensions which again confirms the glassy nature at zero 
velocity in d < 4. The random-force component of the 
disorder correlator does not always vanish in the limit 
w ^ 0, it even diverges for d < 4 in a similar way like the 
KPZ coefficients. 



I. Roughening by disorder 

In the absence of disorder the displacements of the vor- 
tices are isotropic and scale like WxTiY^t) ~ Wyy{r,t) ~ 
w{r,t) ^ b'^''w{r /b,t/b^) with a dynamical exponent 
z = 2 and a thermal roughness exponent C = (2 — d)/2. 

The effect of disorder on the correlations can be esti- 
mated in the most elementary approximation as follows: 
Assuming that the typical displacements are "small" we 
might simply neglect them in the argument of $ in 
Eq. (H). More precisely, the validity of this approxi- 
mation requires Ui — U2 to be small in comparison to 

In this approximation the action is still bilinear in the 
fields and of the same functional form (|l^) as in the pure 
case. From A^"^ originates an additional contribution 



Ar(q,c^) 



^*(k)(5(tj + k-v) 
Q 



(67) 



which is the pinning force correlator as "seen by the per- 
fectly ordered vortex lattice" . Here k is k = Q -I- q. 

The large scale properties of the vortex lattice are in a 
good part governed by the behavior of AF at small q and 
Lu. Note that there is an important difference between the 
properties of components ATxx{(l,^) and ATyy{q,Lu). It 
follows straightforwardly from Eq. ( |67|) that the corre- 
lator $Q^(k) contributes to the asymptotic behavior not 
only at k = but at aU RLV k = Q with Q-v = 0. There- 
fore Arj,j^(q, cj) shows the same qualitative behavior at 
small q and lu for both random potential (similar to un- 
renormalized disorder) and random force case (appearing 
upon coarse-graining). In both cases AF^j^ S{uj + q - v) 
with a finite prefactor for q_,Lu 0. In contrast, the 
behavior of the AF2;a:(q, uj) generated by the random po- 
tential is qualitatively different from that of a random 
force. Namely, in the random potential case ATxx (q, ^) 
vanishes as Arxx{<l,oj) ~ 9x'^(^ + Q • v), whereas for 
the random force case the (5-function again has a finite 
prefactor. Thus the random-force character of the coarse- 
grained disorder will change the asymptotic behavior of 
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those quantities which depend on T^x but will leave in- 
tact those depending on Tyy only. 

In the Gaussian approximation disorder does not mod- 
ify the propagator T and, according to Eq. (|20|), leaves 
also the response function unchanged. However, the ad- 
ditional contribution (|67| ) to the propagator f generates 
an additional contribution 



Xc = rjvyl/c 



(70b) 



^ "'-'^''Kj i^ + kxv) (68) 
c q 



to the correlation function. 

It is interesting to examine the dependence of correla- 
tions on space and time. Since disorder is fixed in the 
laboratory frame, one might expect that AC is also sta- 
tionary in that frame, i.e. ACap{r,t) — Cap{v + vi,0) 
or Cafi{(\,uj) oc 5{uj + qxv). An inspection of Eq. (^) 
immediately reveals that only RLV with ~ give 
such stationary contributions. All other contributions, 
which reflect the discreteness of the vortex lattice in the 
direction of the velocity are not stationary^ neither in the 
laboratory nor in the comoving frame. 

Using the unrenormalized disorder <I>Q/3(k) = 
/Ca/c;3Aoe-5«''=' Eq. (|8|) implies (in d < 3) 

AM^,,(r, 0) « for r » a , (69a) 



that limit the validity range are obtained hy w = W^a — 
(generalizing the static Larkin length) have been in- 
troduced by Giamarchi and Le Doussal |3l| ]. 

If one takes into account that a random-force is 
generated, which we may approximate as $0^(1^) w 

'^1^/3^"^^^'^^ with the coefficients from Eq. (6^) [an expo- 
nential decay for large k follows from Eq. (|)], then also 
the x-component of the displacement becomes rough in 
d < 3: 



AI^,,(r,0) 



AVK,,(r,0) 



a?4)x 



cx 



for rjvy'^ ::§> c\x\ , (71a) 
for rjvy'^ <^ c\x\ . (71b) 



AWyyir,0) 



a^Ao 



3-d 



— for rjvy^ > c\x\ , (69b) 



cx 



for rjvy'^ c\x\ . (69c) 



AWxx has contributions from the vicinity of all RLV to 
order . AWyy is dominated on large scales by contri- 
butions with RLV Qx = 0, whereas RLV with Qx ^ 
give only finite contributions as to AWxx- As for C 
also only the contributions Qx = are stationary in the 
laboratory frame. This includes all contributions that 
roughen the VL on large scales. 

In this perturbative result the transverse displacement 
components exhibit much stronger fluctuations than the 
longitudinal component. An anisotropy emerges requir- 
ing thus a distinct scaling for the displacement com- 
ponents parallel and perpendicular to the velocity and 
also a distinct scaling for their dependence on coordinate 
distances parallel and perpendicular to v. While Wxx 
does not reveal a well-defined scaling behavior at large 
scales since finite temperature and disorder give finite 
contributions (in d > 2), the scaling of Wyy is domi- 
nated by the disorder contribution with Wyy{x,y,t) ^ 
h^-'^Wyy{{x + vt)/b^,y/b,0). The divergence of Wyy on 
large scales indicates that the Gaussian approximation 
loses its validity at large scales, since the initial neglect 
of the dependence of the force correlator on the displace- 
ments breaks down. The characteristic length scales 



a?Ao 



(70a) 



However, the fluctuations of the transverse component 
are eventually more pronounced since (j)yy and "dyy di- 
verge on largest scales (for ~ A*^ 0). This di- 
vergence is logarithmic in d = 3 and algebraic in d < 3. 
Because of this divergence renormalization effects are ex- 
pected to modify the roughness exponents found per- 
turbatively. In addition, one would expect naively that 
AWyy increases on large scales even faster than in Eq. 
(|69|). However, for consistency one should take into ac- 
count not only the one-loop corrections to F, but also to 
r. In the static case, where the Gaussian approximation 
yields a roughness exponent Co = (4 — d)/2, the actual 
roughness is only logarithmic, ^ = 0. There the disor- 
der contribution to the force correlator F, that tends to 
increase the roughness, is balanced by the contributions 
to r, notably an increase to the elastic constants that 
increases the stiffness of the lattice. 

For many disordered systems different approaches to 
treat a balance between several diverging terms have 
been fruitful. One possibility are self-consistent ap- 
proaches like that of Sompolinski and Zippclius |47| ] 
which treats the coupling of modes on a mean-fleld 
level. This approach has been applied for example to 
spin glasses and elastic manifolds |4^. However, 
since typically the self-consistency takes into account 
only first-order corrections A^^'^ to the action, it can pro- 
duce only approximate values of scaling exponents. In 
the present situation these approaches would completely 
miss the physics arising from the divergence of (pyy. It 
would therefore be necessary to extend this scheme to 
the second-order corrections which make a solution of 
the self-consistency equations even more involved. An- 
other approach is provided by the renormalization group 
(RG), that in principle can be extended systematically 
to arbitrary perturbative order. 

A consistent and reliable treatment of the large scale 
properties requires the simultaneous handling of sev- 
eral complications: (i) The anisotropics as discussed 
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above start to interplay under RG interaction with the 
anisotropics of the friction coefhcient, stress coefficients, 
elastic constants, and KPZ terms, (ii) Since disorder 
roughens the VL, as seen already within the Gaussian 
approximation, the generated KPZ terms are more rel- 
evant than the elastic couplings according to scaling ar- 
guments. Therefore, a priori they also need to be taken 
into account and are expected to modify the large scale 
physics qualitatively as soon as the medium is rough. 
The relevance of anisotropics in the KPZ terms in the 
absence of disorder has been studied for single driven 
vortex lines by Ertas and Kardar , who find a variety 
of different physical regimes depending on the anisotropy 
of elastic constants and KPZ coefficients, and by Hwa |5^] 
for driven line liquids. For CDW in d = 1, 2 Chen et al. 

find non-trivial scaling exponents (i.e. differing from 
the scaling found in the Gaussian approximation) due to 
the KPZ terms also in the presence of disorder. 



V. DISLOCATIONS 

Our discussion has been restricted so far to the elastic 
approximation, neglecting topological defects (like dislo- 
cations) in the vortex lattice. Upon increasing velocity 
all effects of disorder become weaker, and we expect also 
the length beyond which dislocations become relevant to 
increase and even to diverge for temperatures below the 
melting temperature of the pure system. Then the in- 
teresting question arises: what kind of defects lead at 
smaller drift velocities to a destruction of the coherence 
of motion, and what kind of order can survive? 

Balents, Marchetti, and Radzihovsky (BMR) ||2| have 
proposed the existence of a smectic phase where vortices 
are correlated over much larger distances parallel to the 
drift velocity than perpendicular to it. 

We propose here a picture for the formation of a smec- 
tic phase within a phenomenological approach based on a 
generalized Lindemann criterion, examining the relative 
fluctuations in the distance of neighboring vortices. In- 
stead of addressing the topological defects explicitly, we 
rather focus on their effect, namely the destruction of the 
neighborhood of vortices. 

Since the positional fluctuations of neighboring vor- 
tices are a small-scale feature, it is sufficient to take into 
account disorder within the Gaussian approximation (i.e. 
neglecting the dependence of the pinning force on the 
vortex displacements) , the large scale properties of which 
have been examined in Eq. ( |69| ) above. The disorder con- 
tribution to the relative displacement of two neighbored 
vortices ri and r2 separated by a bond vector a = ri2 is 
given by 



Au)(a) 



/ (l-e'i-) AC„„(q,c^) 
(a-q)2fc2A(k) 



with AC taken from Eq. (pq). 

Now one can compare the shaking of a bond parallel 
to velocity (a — ax) and "perpendicular" to the veloc- 
ity (a = ay). Strictly speaking, there are no bonds with 
a _L V in a hexagonal lattice, by "perpendicular" we mean 
the bonds making the 60° and/or 120° angles with the 
velocity. This simplified treatment does not change our 
qualitative conclusions altering slightly only the unim- 
portant numerical factors. In the limit of large velocities 
we find 



Aw(ax) 



Aw(ay) 



oAo 



(73a) 
(73b) 



Q 



(72) 



From this result one sees immediately that in this limit 
the bonds perpendicular to v experience much stronger 
shaking effects than the bonds parallel to velocity. Conse- 
quently, these perpendicular bonds linking different chan- 
nels are expected to break more easily than the parallel 
bonds. This implies that the vortex structures has much 
longer correlations in the direction parallel to the veloc- 
ity than in the other directions, in agreement with the 
anisotropy of the dynamic Larkin lengths ( [70| ) . This re- 
sult leads to the phase diagram depicted in Fig. |^ that 
has been discussed already in the introduction. 

Ultimately it is desirable to have a more systematic ap- 
proach to the effects of dislocations in the driven medium. 
A first step in this direction, the study of the dynamics of 
single dislocations, was made in Ref. In order to de- 
cide whether free dislocations destroying the topological 
order of the lattice are present, it is necessary to study 
the dynamic stability of dislocation loops (in d = 3) or 
of dislocation pairs (in d = 2). In particular in d = 2 
one can expect a description of the dynamic phase tran- 
sition in terms of the Kosterlitz-Thouless transition ]5^ ] 
generalized to non-equilibrium systems. 



VI. CONCLUSIONS 

We have constructed a general approach to the driven 
dynamics of dirty periodic media based on the MSR tech- 
nique. The developed scheme provides a regular and con- 
sistent derivation of the effects of disorder on the sliding 
motion. At present, however, we have restricted ourselves 
to the second order dynamical perturbation theory, yet 
sufficient to draw several fundamental conclusions con- 
cerning the high-velocity behavior. 

We have derived the coarse-grained equation of mo- 
tion (|l|) for periodic media in the presence of disorder. 
We have found a renormalization of the bare system pa- 
rameters like friction coefficients, elastic constants, and 
the friction force within the one-loop approximation as 
well as new couplings giving rise to the disorder induced 
stresses, KPZ non-linearities and an effective disorder 
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with a random-force character evolving from the orig- 
inal random potential |30|. The presence of such terms 
has been proposed for one-component systems like driven 
interfaces js^] and CDW without analytic derivation. 

The appearance of divergent parameters under coarse- 
graining is much more subtle in the driven system than in 
the system in equilibrium. For complete coarse-graining 
(A*~ — > 0) we found a divergence in the correlator of the 
transverse components of the effective thermal noise ^* 
and the random pinning force f* in d < 3. Thus in the 
driven state the upper critical dimension is reduced by 
one and there are less divergent parameters in compari- 
son to the static case, where also the friction coefficients 
and elastic constants diverge already in the first order of 
perturbation theory for d < 4. 

The divergences in the correlator of the (persistent and 
non-persistent) random force components perpendicular 
to the drift velocity appear only for systems with a peri- 
odic structure transverse to the velocity. Therefore there 
is a fundamental difference between the dynamic behav- 
ior of CDW, which have only one component, and VL 
with more than one displacement components. 

The standard way to test glassy properties of the sys- 
tems in question is to examine the large scale behavior 
of the disorder-induced corrections to the physical quan- 
tities, the most marking of which is the friction coeffi- 
cient. Its divergence is immediately related to an ex- 
tremely slow dynamics that is dominated by infinitely 
high barriers. A large-scale divergence of the perturba- 
tive corrections would then indicate the glassy behavior. 
However, a divergence of the first order correction to the 
friction coefficient is absent in the driven case. 

Since the vortex system is already in motion, the fric- 
tion coefficient rj*^ , which describes the velocity response 
for a change of the amplitude of the driving force, has 
to be finite since the drifting lattice already overcomes 
the potential barriers in this direction. Nevertheless, one 
could expect that the friction coefficient ri*y, which de- 
scribes the velocity response for a change of the direction 
of the driving force, could still diverge due to infinite bar- 
riers for a transverse motion of the lattice. However, since 
rjyy is finite, these barriers can only be finite. This im- 
plies a linear transverse transport characteristic for small 
transverse forces at finite temperatures. It is still possi- 
ble that a true critical force exists at T = 0. But this is 
not a signature of glassy behavior. An instructive com- 
parison is provided by a single particle in a sinusoidal 
potential, which has no glassy properties jss) . Its trans- 
port characteristics has a finite critical force at T = 0. 
For low temperatures it has an exponentially large but fi- 
nite friction coefficient at small velocities (due to thermal 
activation) that crosses over to the smaller bare friction 
coefficient at large velocities. 

Despite of the absence of generic glassy features like 
the existence of transverse barriers, we find a tendency 
to form channels directly from the presence of the stress 
coefficients (|6^) in the coarse-grained equation of motion. 

To discuss the physical meaning of the effective tem- 



perature 'd* , it seems appropriate to emphasize that it 
is actually the weight of the noise correlator. In a sys- 
tem at equilibrium (and for our bare system), ?? is pro- 
portional to the product of the temperature T and the 
friction coefficient 77. Since we have found a finite correc- 
tion to the friction coefficient, the divergence of can be 
interpreted as a divergence of the effective temperature. 
However, in the non-equilibrium situation under consid- 
eration, there is no well-defined meaning of a "tempera- 
ture." One can speak only about an (non-unique) effec- 
tive temperature if one specifies what physical property 
of the non-equilibrium system is compared to an equi- 
librium system. Since the divergence of t?* arises from 
fluctuations on the largest length scales, only the degrees 
of freedom on asymptotically large scale can be related 
to an infinite effective temperature. This means that the 
driven lattice in the presence disorder is on large length 
scales much more rough than the lattice in the absence 
of disorder. 

In general, divergent parameters signal a break-down 
of perturbation theory at large scales. Therefore the 
question about the asymptotic large-scale behavior of the 
system can be conclusively addressed only by a system- 
atic RG treatment that includes implicitly all orders of 
perturbation theory, which goes beyond the scope of this 
work. 

In the absence of a formal derivation, we propose the 
following speculations regarding the existence of (quasi-) 
long-range order at highest driving forces. It is essen- 
tial to distinguish between CDW-like systems (having 
only one "displacement" component) and VL-like sys- 
tems (having more than one "displacement" component) 
because of fluctuations in the displacement components 
perpendicular to the velocity. In the CDW case, where 
there is only the component parallel to the velocity, it 
has been argued in Ref. ||3^ that the random forces 
along that direction lead to a roughness with an expo- 
nent C = (3 — c?)/2 that is not reduced by renormalization 
effects on the largest scales. 

We believe that the situation could be different in 
the case of VL. We have shown within the perturbative 
framework that the perpendicular displacement compo- 
nents fluctuate much stronger than those of the paral- 
lel component and are subject to a diverging random 
force correlator. In this case one has to take into ac- 
count that the strong perpendicular displacements wash 
out not only the perpendicular components of the pin- 
ning forces but also the parallel components and there- 
fore qualitatively reduce the true large-scale roughness in 
all directions. 

This speculation is formally supported by the struc- 
ture of the perturbative corrections obtained in Section 
III. To be specific, we discuss the random force correlator 
Eq. (^). It contains exponential factors exp{— k-W-kj 



that have been neglected in the evaluation in Section IV 
This is legitimate in lowest order perturbation theory, 
where only thermal fluctuations contribute to W. How- 
ever, under iterating the perturbative expressions (this 
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is essentially the idea of an RG), one should take into 
account also the disorder contributions to W. Since we 
have found that disorder roughens the system, W(r, t) di- 
verges on large scales and the exponential factors decay 
for large k, r or t, suppressing the corrections in com- 
parison to the lowest order perturbative results. In this 
way all corrections generated by disorder are suppressed, 
except from the first order correction (|4l| ) that tends to 
compensate the disorder itself. 

One also immediately recognizes that the roughness of 
the components perpendicular to the velocity (the large 
scale behavior of Wyy) influence e.g. the random force 
correlator cjjxx in directions parallel to the velocity and 
vice versa. 

This mechanism is the same in all perturbative expres- 
sions and, as we believe, coud reduce the true roughness 
qualitatively in comparison to the perturbative result in 
2 < d < 3 for the largest velocities. 

In addition to the aspects of the high-velocity phase 
we have evaluated our general perturbative results also 
in the limit of vanishing velocities, when the depinning 
transition is approached. We found disorder to be rele- 
vant in c? < 4 in agreement with Larkin's original static 
analysis Q and dynamic approaches to depinning as a 
critical phenomenon [p||2^]. As an additional feature 
our perturbative analysis revealed the relevance of the 
non-equilibrium contributions (diverging KPZ terms and 
random-force correlator) to the equation of motion. 

To be specific, let us consider the example of the KPZ 
couplings as true non-equilibrium couplings. It seems 
surprising that they do not always vanish in the limit 
w — > 0, where one naively expects the FDT to hold and 
all non-equilibrium terms to vanish. To resolve this para- 
dox note that the zero- velocity limit has to be taken with 
care since it does not commute with the limit —> 0. 
The observed divergence occurs only if one takes A^ 
before w — > 0, since the divergence arises from the in- 
frared contributions to the integration over q. Physically 
these contributions are related to the diverging energy 
barriers on large scales. These diverging barriers imply a 
diverging relaxation time and persistent memory effects 
of the system, which are the origin of the survival of non- 
equilibrium terms. In other words, in a glassy system 
relaxing from a (globally drifting) non-equilibrium state 
into its equilibrium state after switching off the current, 
there will be regions which are still drifting, their dy- 
namics being governed by an effective non-equilibrium 
equation of motion. If, on the other hand, one considers 
the KPZ couplings for finite A^, they are finite and van- 
ish ~ V for velocities v ^ c^A^^/yy. For this reason these 
terms have not been taken into account in the previous 
studies of the depinning transition. However, the ober- 
vation that these non-equilibrium terms diverge at small 
velocities for A< = could indicate that even for weak 
disorder the depinning transition is rounded in a very 
narrow region by rare plastic effects |5^ , which are not 
captured by the phenomenological Lindemann criterion. 



In the final stage of preparing this manuscript we be- 
came aware of a preprint by BMR that extends their 
earlier reference jS^ and where the authors come to sim- 
ilar conclusions. In addition to perturbative results a 
RG analysis of a simplified "toy model" for the trans- 
verse displacement component has been performed, but 
it does not seem to yield results that differ qualitatively 
from those obtained by the perturbation theory. 
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APPENDIX A: NOTATION 

For the convenience of the reader we summarize here 
our notation. The total dimension is d = dy + d± with 
the inner vortex dimension and the number d± of 
displacement components. Latin indices run over all d 
components, whereas Greek indices run over d± com- 
ponents. We use summation convention for indices, i.e. 
k • k = fcafca = k^. 

Some short notation: 

ri2 ri - r2 , ti2 ti - t2 

Si2 ■■= <5(ri2)<5(ii2) , <^(r) a'^^S^Jir^^) 



Ji Jvi Jti Jt J Jr j.^ J 

d'^q 



'ti 

d'^k 
(2^ 



IBZ 



duj 
2^ 



For discrete space-components the integration has to be 
replaced by a sum with a factor a'^^ , the "volume" per 
vortex. In dj^ = 2 one has — ^o/B. q-integrals run 
only over the first Brillouin zone (IBZ) but k runs over 
whole momentum space. 

For two-point quantities a conjugation is defined by 

Gt,/3(k,c^) :=G;3a(-k,-c^) , (A2) 
ij = G^c(^(ri — Tj, ti — tj) :— Gpai^j — fi, tj — ti) . 
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APPENDIX B: SECOND ORDER CORRECTION 
TO THE DISORDER CORRELATOR 



In order to find the second-order contributions to the 
disorder correlator, one has to look at terms of order 
O(w^) in (|2^), which we separate into two contributions 
(all superscripts are dropped again) 



^g^=^r+-^r (Bl) 
specified below. A further condensation of the products 

— e*-^i^ (B2) 

and <& 12Q/3 := ^21/3q helps us to keep the structure of 
the expressions transparent: 



Ag''^ = / ([ill . $5^; . (Gtgi . k' + (G+32 - Gt34) • k")][(k' • (Gi4 - G34) + k" • G24) • ■ O12] e-^ 

^ Jl234k'k" I 

-[ui . • Gtai . k'][k" . G24 ■ *^4t" ■ fi2] , (3a) 

4''"^ = -7 / ["1 • *i2 • "2](k' • (Gi3 - G23) - k" • G43) • • ((Gt4i - Gt42) • k' + Gt43 • k") e"^ (3b) 

4 Jl234k'k" 

with the exponentials (generalize Debye- Waller factors) 

A EE A^V' 1 ^[k' • (Ui - U3) + k" • (U2 - U4)]2) 

= k' • Wi3 ■ k' + k" ■ W24 • k" - k' • W12 ■ k" + k' • Wi4 • k" + k' • W32 ■ k" - k' ■ W34 • k" , (4a) 
B EE B^;^' := A5^V' . (4b) 

The second term in the curly brackets of Eq. (^) arises from the subtraction in the definition of the cumulant 
(ylP'" ylP'")c := (^P'"^P™) — (^P'")^, which correspond to "disconnected diagrams" in field-theoretical language. In 
Eq. (||b) one further term ~ G^34G34 actually vanishes due to causality, i.e. G''^34 ~ 6(^4 > t^) and G34 ~ 6(^3 > t^). 

Contributions that represent effective disorder are identified as follows: in (||) the response fields are evaluated at 
points {ri,ti) with i = 1,2. Disorder contributions are those which persist for |ii2| — *■ 00. This requires that Ui and 
U2 are not connected (even indirectly) by response functions. In this limit the factors and simplify, since 
Wjj Cq if the points i and j are unconnected. Therefore we are left with 



2 



1234k'k" 

-W34)-k" 



X {(Gtgi . k' - Gt34 . k")(k' • Gi4 - k' • G34) e-k" Co k"-k'.Wi3 k'-k'.(Wi4- 
+ (Gt32 • k" - Gt34 . k")(k" • G24 - k' • G34) e-k'-Co-k'-k".W24-k"-k".(w.3-W43)-k' 
+Gt32 • k"k' • Gi4 e-('''-''")-c-('''-''")-'''-(w-+w-)-''" } . ^t;^''" . U2 . (5a) 

Terms ~ Gt3iG24 are non-persistent and therefore do no longer appear in Vanishing terms ~ G^34G34 have 

been inserted by hand to complete the squares. The second contribution becomes analogously 



4 



1234k'k" 



({(k' • Gi3 - k" • G43) • *3r • (Gt4i . k' + Gt43 . k") e-k'-Co-k'-k".W34-k"+k'.(Wi3-Wi4)-k" 

-(k' . G23 + k" • G43) • *3r • (0^42 • k' - Gt43 . k") e-k'-Co-k'-k".W34-k"-k'.(W,3-W,4)-k" 

-2k' . G23 • • Gt4i • k' e-(k'-k")-Co-(k'-k")-k'.(w.3+Wi4)-k" I _ (5b) 
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Here the number of terms was reduced using the relabel- 
ing symmetry 3 <-> 4 with k" — k". 

To identify the corrections to the disorder correlator 
we have to analyze the dependence of the functional on 
the displacement field and the response field. Since we 
are interested in the large-scale physics, we can consider 
the response-functions as local in space and time. The 
labels have been chosen such that u is attached to 1 and 
2. Since in all expressions two response- functions are in- 
volved that connect point 3 and point 4 to point 1 or 
point 2, U3 and U4 may be replaced by the correspond- 
ing Ui or U2. These replacements can be considered as 
lowest order of the expansion (^9|) . 

To give an example we discuss the first term in the 



curly brackets of Eq. (^a). There the response func- 
tion connect points 3 and 4 with 1, whereas 2 is free. 
In this case it is convenient to replace U3 and U4 by ui. 
The exponential factors implicit in the factors then 
can be rewritten as e*'-('"i3+'^*i3)+»k"-(r24+vt24+u2i) ^ 



-jk"-(ri2+vti2+ui2)+jk'-(ri3+vti3)+jk"-(ri4+vti4) 



Now the 



integrand has the same functional dependence on the 
fields as the original disorder action and a contri- 
bution to the correlator correction can be identified. In 
the same way one can proceed with the other terms of 



Eq. (h) and finds (abbreviating R! 



(0) ._ 

ij 



*(2)(k) = -e 



-k Cn k 



34k'k" 



c{*(k') • [(Gtgi • k' - Gt34 • k")(k' • Gi4 - k' • G34) e''''-Wi3-k'-k'.(Wi4-W34)-k" gik'.R5?)+,k".R<«>^(j^ ^ j^//) 

-(Gt32 • k" - Gt34 . k")(k" . G24 - k' • G34) e-''"-^--''"-''"-^^--^-)-''' e'k'.<3'+»k".R<°'^(]^ _ 

-G^32 ■ k"k' • Gi4 e-k'-(Wi4+W32)-k" g»k'.R<«)+^k".R(°)^(^ _ j^, ^ . ^t(_k") 

*(k') [(k' • Gi3 - k" • G43) • *(k") • (Gt4i . k' -t- Gt43 • k") e-k"-W34-k"+k'.(Wi3-Wi4)-k" e»k"-RS' j(k - k') 
-(k' • G23 + k" • G43) • *(k") • (Gt42 • k' - Gt43 . k") e-k"-W34-k"-k'.(W23-W24)-k" ^^"-^'^2 s(^y, _ k') 
-2k' • G23 • *(k") • Gt4i . k' e-k'-(W23+Wi4)-k" g*".(R(°)-R(°')^(]^ _ j^, ^ I 



which in fact does not depend on points 1 or 2, which 
can be eliminated by substitutions for the points 3 and 
4. 

The symmetry $(k) = €>^(k) given for the original 
disorder correlator is preserved after the inclusion of the 
corrections. 

In order to go beyond the locality approximation used 
above, one could include the derivatives of Eq. (p9|). This 
is not done here, since a scaling analysis shows that the 
resulting terms will be less relevant than the disorder cor- 
relator. However, along these lines one can straightfor- 
wardly derive a random KPZ nonlinearity as postulated 
by Krug pi. 
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Y,Z 



l/(3-d) 




l+2/(3-d) 



FIG. 1. Shape of the Larkin domain in the drifting struc- 
ture. The domain is much longer in directions parralel (x) 
than perpendicular (y, z) to v. Its anisotropy increases with 
increasing velocity, see Eq. (|7C|). 
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FIG. 2. Topology of the generic phase diagram of driven 
periodic media (a) for strong disorder and (b) for weak dis- 
order. At largest velocities the medium moves coherently. At 
the solid line a melting occurs into an incoherent phase with 
massive plasticity. This phase can display smectic proper- 
ties and decay at even smaller velocities (dotted line) into a 
fluid-like phase. At even smaller velocities (dashed line) the 
creep phase with metastable states is reached. 



FIG. 3. Geometry of the vortex lattice in d = 1 -(- 2. Vor- 
tices carry a fixed label rj_ and their fluctuating position is 
R. The a;-axis is chosen parallel to the average velocity v. 




V 



FIG. 4. Schematic representation of a moving vortex con- 
figuration. The distance between the solid wiggly lino (posi- 
tion of vortex in perfect lattice) and the dotted straight line 
(actual vortex position) represents the component Uy of the 
displacements field. The force dynamically generated force 
fv ^ XyyxdxUy, visualized by empty arrows at two represen- 
tative positions, makes a vortices acquire the displacement of 
the preceding vortex. 
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